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^~~> . We improve 0(p*) calculation in U{2)l x f/(2)_R chiral theory of mesons by resummation calcu- 

' lation for vector mesons physics and restudy decay p — > tt-k. A complete and compact expression for 

y—i , fpinrip^) (up to 0(p°°)) is obtained, from which an important non-perturbative conclusion is given 
based on convergence and unitarity consideration. 
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(N : 

' The Nambu-Jona-Lasinio (NJL) model [1] and its extensions are widely used to understand hadron physics (for 
\ recent reviews see e.g. [2]). There are various methods to parametrize the NJL version models, among which is the 

■ U{Nf)L X U{Nf)ji chiral theory of mesons [3] with Nf = 2 or 3 (proposed by Li and called hereafter Li's model). This 
I model can be regarded as a realization of chiral symmetry, current algebra and vector meson dominance (VMD). It 
, provides a unified description of pseudoscalar, vector, and axial-vector mesons, which are introduced as bound states 

of quark fields. The basic inputs of it are the cutoff A (or g in [3]) and constituent quark mass m (related to quark 
condensate). This theory has been studied extensively [4-9], in particular it has been used recently [10] to analyze 
the data of g-2 of muon reported by CMD-2 group [11]. Hence Li's model is a good phenomenological model. 
Qh! So far Li's model is only investigated up to O(p^) [3] in perturbative manner. This treatment is reasonable for 
^ ■ chiral perturbation theory (ChPT) [12], because the typical energy there is much less than the energy scale of chiral 
symmetry spontaneously breaking (CSSB) Acssb ~ 27rF^ ^ IGeV. This is not the case here because the typical 
energy p of vector mesons is comparable with Acssb- For example, at p energy scale, p ~ and the perturbation 
parameter p/ Acssb ~ 0.7. Therefore, studies in ref. [3] should be improved by going beyond 0{p'^). But we know that, 
in usual perturbative treatment, it is difficult to deduce effective meson action at higher order. This is a paradox. As 
' a solution to it, we develop a method to obtain the sum of all terms of the p-expansion [13]. We call it resummation 
study. In this paper, we want to use it to restudy typical decay p — > vrvr in Li's model. Other processes can be 
restudied in a similar way. 

Before this study, let us outline procedures of perturbative calculation and our resummation calculation. Given 
a quark model which describes interactions of constituent quarks q and mesons M and satisfies chiral symmetry 
requirement, we can write the Lagrangian as £ = C{q,q,M) = qVq, where V = 'D{M) is characteristic of the model. 
The first step of perturbative calculation is using the method of path integral to integrate out the quark fields: 

e-Jd-'^^MM)^ f[dq][dqyl'''''^^''. (1) 
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After ftmctional integration, the effective action Scs{M) of mesons is formally obtained as 5off (M) = In detD(M). 
To regularize this determinant, Schwinger's proper time method [14] or heat kernel method [15] should be employed. 
This directly results in an expansion in p, and usually up to 0{p^). Attempting to calculate terms of higher orders 
will encounter great difHculty. 

On the other hand, to perform resummation study, we calculate effective meson action via loop effects of constituent 
quarks. This is equivalent to integrating out quarks in path integral, but via it we can obtain effects from all orders of 
p-expansion. Specifically, we divide £ into free- field part £o and interaction part £i, and turn to interaction picture. 
The effective action can be obtained as 

e^s=« = <0|r,e*/'''"'^'|0> 

= / c^V^;^---^;^^^(Pi+P2 + ---+Pn)n„(pi,---,p„), (2) 

where Tg is time-order product of constituent quark fields, n„(pi, • • • is one-loop effects of constituent quarks 
with n external fields, Pi,P2r " ^Pn are their four- momentum. Getting rid of all disconnected diagrams, we have 

SeS = E ^" = / "^"P^T^ ■ ■ ■ T^^^iP^ +P2 + -- ■+Pn)K{Pl, ■ ■ ■ ,Pn), (3) 

where c denotes connected part. Obviously, in eq. (3), the effective action ^eff is expanded in number of external vertex 
and expressed as integral over external momentum. Hereafter we shall call this method proper vertex expansion, and 

call Sn n-point effective action. In terms of proper vertex expansion, the effective actions include informations from 
all orders of chiral expansion. That is, we can do resummation of momentum expansion by this method. That is 
what we need. 

It is instructive to compare these two kinds of calculations. The perturbative calculation can deal with all kinds of 
meson interactions at low orders, while resummation calculation deals with specific process but can include informa- 
tions of all orders. In the sense of p-expansion, resummation calculation is non-perturbative. As we shall see later, 
this character will give important result. 

Having showed outlines of resummation calculation, we turn to Li's model. Because we only focus on p physics, 
the Nf = 2 part of Li's model, i.e., U{2)l x U{2)ii chiral theory of mesons [3] suffices to study decay p tttt. This 
Model is constructed through U{2)l x U{2)r chiral symmetry and minimum coupling principle, and the ingredients 
of it are quarks {u and d), pseudoscalar mesons (tt and rj {u and d component)), vector mesons {p and w), axial-vector 
mesons {a\ and /i(1285)), lepton, photon and W bosons. For our purpose, we just write the relevant Lagrangian: 

£ = q{ip + f +475 - muix))q + ^ml{< V^V^' > + < A^A^' >), (4) 

where = t^o^j + /^j, = Tip^^ + ujf^, u = e'*^"^ , <I> = TiTr' + 77, and < • • • > denotes trace in flavor space. The quark 
part of this Lagrangian can be divided into two parts. The free-field part is £q = q{ip — m)q, and the interaction part 
is 

£l=q{f+4l5-im<i>j5 + \rn^^)q, (5) 

where terms involving more <l>s has been omitted, because they have nothing to do with decay p —> nn. All subsequent 
calculations are performed at chiral limit. 




(a) (b) 

FIG. 1. Calculation of < CD > terms (C, D — V, A, $) from quark loops, a) Two-point 
diagram of quark loop, b) One-point diagram of quark loop for < > term. 
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Before studying decay p — > tttt, we have to do something to obtain physical fields. We should calculate effective 
actions for < VV >, < AA >, < >, < VA >, < > and < A$ > through integration of two-point quark loops 
like Fig. la (where C,D = V, A, $) and one-point quark loop Fig. lb. These actions are calculated as 

2 J (27r27r)4 

( < V^{p)V4-p) > TT[j''SF{k-p)YSF{k)]+ < A^{p)A^{-p) > TririsSpik - p)r75SF{k)] 

- < ^pM-p) > {m^Tr[j5SF{k-p)j5SF{k)]+mTi[SF{k)]) 

~2im < V^{pM-p) > TT[j^SF{k-p)r,SFm - 2im < A^{p)^-p) > Trlji'jsSFik - phsSF{k)] 

+2 < V^{p)A,{-p) > Tv[j''SF{k-p)ri5SF{k)]), (6) 

where, for < > term, we have added contribution from one-point diagram (Fig. lb) of quark loop. This is so 
because is realized nonlinearly. After integrating quark loops, we get expressions including all terms up to 0(p°°). 
But we are not interested in high order terms now. Discarding them and adding the non-quark part of eq. (4), we 
obtain 



1 



_d_^p 

(2 



^ (< V^{p)V^{-p) > {g-'ip^p^ - g^'^p') + ml) 



+ <Mp)M-P)>[-9'{i--^)9^r 



6'ir^g'^ 



+ < $(p)$(-p) > -g^m^p^+ < A^{p)^{-p) > &ig^w?p''^, 
where the logarithmic divergence is absorbed by 



Nr. 



(4,^)^/2 



e/2 



r(2-- 



(d: 



(7) 



(8) 



and the quadratic divergence in < > term has been reduced to logarithmic divergence by the identity r(l — d/2) + 
r(2 — d/2) = — 1 -t- 0(e). The first and second terms of eq. (7) indicate needs for rescalings of and (We see 
that axial-vector yl^ does not satisfy gauge invariance. Because dfj_A^ = when A^, is on-shcll. wc shall ignore the 
p^p^ term when redefining A^.) The non- vanishing < A$ > term indicates that there is a mixing between A^ and 
9^$, which should be erased by shift of A^. Therefore, the redefinitions of and A^ are 



A„ 



-1. 



g^l - N,/67r^g^ g 



and we can obtain in passing the mass formula for them: 



rriv = 



,9^ 



(1- 



67r^5 



2 2)"^^ ~ +mv- 



(9) 



(10) 



After redefining and cancelling the mixing, and then making the kinetic term of $ in standard form, we obtain 



1 



6m^ + m?/ ' 



-C'mi + - 1 



2c 



P2 
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or 



c = 



f2 

2gmy ' 



6(1--) g'm' = Fl, 



(11) 



and rescaling for $ 2$/F$. Eq.s (9)-(ll) are exactly those in [3]. 
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FIG. 2. Quark loops for vertex p tttt. a) Two-point diagram, b) Three-point diagram. 

Now we can study decay p tttt. For vertex p — tttt, we need to calculate two Feynman diagrams of quark loop 
(see Fig. 2). The relevant Lagrangians (after redefinitions) for the first and second ones are 



= ^i-J + |J *^)., = ,( V - ^^7s - f^$7s)a (12) 



respectively. For the first one. c[^^ gives integration of the quark loop as J d'^kTilj^Sp^k — p)SF{k)] oc p^^ (p is 
momentum of vector field V"). Because d^V^ = when is on-shell, contribution from the first diagram vanishes: 

^2 = 0. For the second diagram, the Lagrangian C\ determines that the effective action for this three-point diagram 
is calculated as 



4i7Vc f d^pd^qid'^q2 



Ne f 

p2 / 
(J) J 



gFl J (27r27r)4 



^'^{P + qi+ q2) < V^{p)^qi)^{q2) > 



[^qi.q2pTr[^''SFik + qi)Yl5SF{k)Yl5SF{k - 92)] + m^Tr[^''SF{k + qxhr,SF{k)^nSF{k - 92)] 



(27r)4 \g 

-—q^,T^[YSF{k + qih^i^SFikh^SFik - q2)] - —q2.TT[j^SFik + qih^SFik^'l^SFik - 92)]) ■ (13) 

After integrations of quark loops and considerations of d/^V^ = and chiral limit 9i = 52 = = 0) the effective 
action for vertex p — ww becomes 

Sp.. = I ^^A^^^^p»^'i-p - Q^'m-iin/p^Ap^ (14) 

where 

Mp^) = [m2(l8(l - 2c/g)n^g^ - (2 - 3c/g)N,) -p\c/gf{Qn^g^ - N,) 

/ xdx / dy[^ — -\'m?{l-2c/g + xy) 

ttWwo Jo ym? -p^x{l-x){l-y)V 

+p'^(x{l -x){l- y){l + xy) - c/g{l + 2x- 2x^ - 3xy + 2x'^y) + c^/g'^{l - xy)^ 

2 

-ln(f - ^x{l - x){l - y)j (m^il - Ac/g + 3xy) - cV/5^(l - a^y)) 
_ 12(iVe + Sg^TT^) - (24c/g)(2iVe + 3g^7r^) + 40N^{c/g)^ 2 c^(10A^c + 36g^7r^) 2 



4iVe(3-I2c/5+10(c/5)2)m2-47Ve(c/5)V Um^ - p\ , / p^ 

'V4^^^" ^ ^ 

The factor fp-rnrip^) is complicated, expansion of it in p up to 0{p^) (corresponding to expansion of SpT^v up to 
0(p8)) is 



2 
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N,{1 - 2clgf - 12^TT^ 2 , (1 - 4c/5)iVc 4 , (l " 4c/5 + {c/gf)N, 



Mp') = - 1 + ^2 + P' + .ool^J;;^ P' + , (16) 



67r2/^ 607r2m2/2 ^ A20'K'^m^f^ 



where condition (II) has been used. As we can see, the sum of the first two terms in this series is just /p^^ in ref. [3] 
when A^c = 3 and p^ = rn^ is set. 

Moreover, expression (15) shows that series (16) is in fact an expansion in dimensionless quantity p"^ = p^/(4to^): 
fp-mrip^) = 2/g -\- c\p^ + C2P^ + c^p^ + ■ ■ ■■ Convergence of this series means < 1 or < 4m^, which at p energy 
scale means the constituent quark mass m > mp/2 ~ 385MeV. The same conclusion can also be obtained by unitarity 
consideration. Unitarity of the theory demands that, at leading order of I/A"c expansion, there should be no imaginary 
parts in transition amplitude of meson decay [16]. The expression (15) for /p^^r is just at the leading order of 1/A^c 
expansion, therefore, it must be real. Thus p'^ < Aw? must be ensured in eq. (15), which means m > mp/2 at p 
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energy scale. Wc should point out that this conclusion is obtained from non-perturbative expression (15) which is 
characteristic of resummation study. 

The difference between our result of m > 385MeV and the one of to = SOOMeV in ref. [3] is understandable. The 
parametrization m = SOOMeV in ref. [3] is consistent with its 0{p'^) calculation, and the phenomenology studies are 
good there. Now that we have resummed all terms in ;3-expansion, using convergence or unitarity analysis, we have 
that TO > 385MeV. In fact, such high constituent quark mass is also reasonable from another aspect. When fitting the 
hadron spectra, we have to adopt high constituent quark mass, or we fail to account for the observed masses of scalar 
nonet [17]. It is argued in ref. [17] that high constituent quark mass might be a consequence of lack of confinement of 
the model and could be avoided if wc knew how to add confining forces to the model. This is suggested by the results 
of constituent quark models which use confining interactions [18]. 

To conclude, we improve perturbative calculation in Li's model by resummation study. We first illustrate proper 
vertex expansion method in chiral quark model, and then use it to perform resummation study on decay p tttt. The 
complete and compact expression of fp-wTrip^) (up to 0(p°°)) has been derived, from which the explicit expression of 
effective action 5^^^ up to 0{p^) is easily obtained. We have seen that this method of resummation study is a powerful 
method to catch informations from all orders of p-expansion, from which we obtain non-perturbative conclusion that 

TO > TOp/2. 
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